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, We describe a model of spectral energy distribution in supercritical accretion disks (SCAD) 

based on the conception by Shakura and Sunyaev. We apply this model to five ultra-luminous 
Oh" X-ray sources (ULXs). In this approach, the disk becomes thick at distances to the center 



less than the spherization radius, and the temperature dependence is T oc r -1 / 2 . In this 
region the disk luminosity is Lboi ~ -^Edd hi(Af/MEdd), and strong wind arises forming a 
wind funnel above the disk. Outside the spherization radius, the disk is thin and its total 
luminosity is Eddington, ^Edd- The thin disk heats the wind from below. From the inner side 
of the funnel the wind is heated by the supercritical disk. In this paper we do not consider 
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Q-H, Comptonization in the inner hot winds which must cover the deep supercritical disk regions. 

Our model is technically similar to the DISKIR model of Gierlinski et al. The models differ in 
disk type (standard — supercritical) and irradiation (disk — wind). We propose to distinguish 

a ■ 

between these two models in the X-ray region ~ 0.3-1 keV, where the SCAD model has a flat 
vF v spectrum, and the DISKIR model never has a flat part, as it is based on the standard 

> 

a-disk. An important difference between the models can be found in their resulting black 

m ■ 

\q ' hole masses. In application to the ULX spectra, the DISKIR model yields black hole masses 



O ■ ~ 10 Mo- 

rn 



of a few hundred solar masses, whereas the SCAD model produces stellar-mass black holes 



1. INTRODUCTION ulation III stars [2] can form binary systems due 

to tidal captures of single stars. However, the 

Ultra-luminous X-ray sources (ULXs) are X- expected frequencies of such t M BHs are not high 

ray sources with luminosities exceeding the Ed- and do not agree ^ the observed ULX f req uen- 

dington limit for a typical stellar-mass black deg [3] _ Mogt of thfi ULXs ^ associated with 

hole [1]; their luminosities are 10 39 -10 41 ergs" 1 . staj formation regions [4] and young star clus . 

The most popular models for ULXs involve ei- terg [5] They are not distributed throughout the 

ther intermediate-mass black holes (IMBH, 10 3 - galaxy ag WQuld be expected for IMBHs originat _ 

10 4 M ) or stellar-mass black holes (~ 10 M Q ) ing from low . me tallicity Population III stars, 
accreting at highly super-Eddington rates. Both 

scenarios require massive donors in a close orbit. IMBHs may be produced in the process of 

IMBHs originating from low-metallicity Pop- runaway merging in the cores of young stellar 
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clusters [6]; in this case they would stay within 
the clusters. Poutanen et al. [5] have found that 
the majority of the brightest X-ray sources in 
the Antennae galaxies are located close to very 
young stellar clusters, but not within them. They 
concluded that the sources have been ejected in 
the process of star cluster formation due to dy- 
namical few-body encounters and that the ULXs 
are massive X-ray binaries with the progenitor 
masses larger than ~ 50 Mq. 

The model where the ULXs are supercriti- 
cal accretion disks (SCADs) requires geometrical 
collimation of radiation into the observer's line of 
sight [7, 8]; this model explains the objects up to 
possibly 10 41 ergs" 1 . It was suggested the ULXs 
are face-on copies of SS433, the only known su- 
percritical microquasar in the Galaxy [9]. The 
binary SS 433 contains a stellar-mass black hole 
accreting in a supercritical regime at ~ 300-500 
Eddington mass accretion rates. In this system 
the wind forms a wide funnel along the disk axis, 
which collimates the disk radiation. In our recent 
optical spectroscopy of ULX counterparts [10], 
we have discovered that their optical spectra are 
very similar to that of SS433. 

If the ULXs are SCADs with stellar-mass 
black holes, one may expect to detect discrete 
lines or edges in their spectra because of the mas- 
sive outflows in the supercritical regime. The X- 
ray spectra of ULXs do not show any emission or 
absorption features; the best quality data [e.g., 
11] yield an upper limit on the equivalent widths 
of a few tens of eV. A common property of the 



ULX spectra is a high-energy curvature [12-14] 
with a downturn observed between ~ 4 and 
~ 7 keV. The curvature suggests that the ULX 
accretion disks are not standard. The inner 
parts of the disks may be covered with an out- 
flow or with a relatively cool and optically thick 
corona [12, 13], which Comptonize the inner disk 
photons. 

The ULXs are frequently surrounded by nebu- 
lae [15, 16]; the nebula shapes are similar to those 
powered by jets or weakly collimated winds. In 
some ULXs, radial velocity perturbations in the 
associated nebulae were detected [17]. 

Strong radiative beaming is not necessary for 
the ionization of the ULX nebulae. Intermedi- 
ate geometrical beaming (B ~ 3-5) is enough to 
produce both the observed luminosities and the 
observed ionized nebula shapes. Besides that, 
not only X-rays may ionize the nebulae. In the 
optical study of MF 16, the nebula surround- 
ing NGC6946 ULX-1, Abolmasov et al. [18] 
have found that UV and far UV radiation is 
able to produce the MF 16 spectrum; they sug- 
gested that the ULX systems may also be ultra- 
luminous in UV. It has been found that at least 
some of the ULXs are indeed very bright UV 
sources [19, 20]. However, both the IMBHs and 
SCADs are expected to be bright UV sources. 

Several distinctive properties of ULXs may be 
explained in these two competitive models. In 
spite of much effort in the ULX studies, the main 
questions still remain open— what are the black 
hole masses in ULXs, do they have standard ac- 
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cretion disks, do they constitute a homogeneous 
class of objects. 

One may suggest that the common proper- 
ties of ULXs 3X6! cl hot disk wind which Comp- 
tonizes the inner disk radiation, making a Comp- 
ton tail in the X-ray spectra, and a very luminous 
UV/optical source. The ULX optical counter- 
parts show hot blue spectra; the absolute magni- 
tudes of the nearest and well-studied ULX star- 
like counterparts (not clusters) in the HST im- 
ages are My ~ — 6 ± 1 mag [21]. 

Both the high UV/optical luminosities and 
the Compton tails may be reconciled with a self- 
irradiating accretion disk. Such a model has been 
suggested by [22, 23], it was originally developed 
for LMXBs. If the disk is flat, as is the case in 
the standard disk model, one finds the temper- 
ature dependence T(r) oc r -3 / 4 , the same as in 
the disks without self-heating. The authors sug- 
gested that the disk, for some reason, is warped. 
In this case, one may expect a flatter dependence, 
T(r) oc r -1 / 2 , providing strong UV and optical 
fluxes from the outer parts of the disk. This 
DISKIR model has been successfully applied in 
several ULXs [20, 24, 25]. 

The SCADs may potentially explain both the 
high UV/optical luminosities and the power-law 
Compton tails in ULXs. The disk wind photo- 
sphere may be a luminous source with a hot black 
body-like spectrum [8, 9]. Deep inside the super- 
critical disk funnel, where the jet is formed and 
collimated, hot winds make favorable conditions 
for the Comptonization of the surrounding pho- 



tons radiated by the funnel walls. The hard radi- 
ation is expected to be mildly collimated by the 
funnel. The outer walls of the funnel constitute 
the SCAD's wind photosphere. The jets may 
power and shape the nebulae associated with the 
ULXs as it is in SS 433. 

The supercritical accretion regime was first 
described by Shakura and Sunyaev (1973) [26]. 
They introduced a "spherization radius" in the 
disk, R sp , below which the disk becomes super- 
critical. The supercritical disk is thick with a 
strong mass loss. In the supercritical region, 
wind inevitably appears; its velocity is close to 
virial, as in stellar winds. At highly supercriti- 
cal accretion rates (Mo >> MEdd); the extended 
wind photosphere (R p h » Rsp) conceals the 
wind formation region, as is observed in SS 433. 

The recent 2D RHD simulations [e.g., 27, 28], 
which take into account both the heat advection 
and photon trapping, confirm the main ideas of 
the Shakura-Sunyaev SCAD model. Poutanen et 
al. [8] have considered ULXs within the frame- 
work of the Shakura-Sunyaev approach; they 
confirmed that the main features of the ULX 
spectra may be formed in the extended photo- 
spheres. 

In this study, we simulate the spectral en- 
ergy distributions (SEDs) of SCADs based on 
the Shakura-Sunyaev approach. We find that 
this model is similar in the main features to the 
DISKIR model [22, 23], however, the models dif- 
fer in many details. In our model, the inner parts 
of the supercritical disk (R < R sp ) heat the outer 
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wind, whereas in the Gierlinski et al. model the 
inner regions of the accretion disk heat the outer 
parts of the disk. In this paper, we have not yet 
included the Comptonization in the model, and 
consider the model to be valid up to < 1.5 keV. 
We plan to introduce the Comptonization in the 
next paper. 

It is important to have UV data points as 
close to the X-ray range as possible, because the 
spectral shape in this region may be critical for 
distinguishing between the models. Photomet- 
ric data in < 2000 A, bands are available only 
for two ULX counterparts. It is also important 
to have simultaneous data. When the variabil- 
ity of the ULX counterparts is not very high in 
the optical range (about 10-20%, [21]), the X- 
ray variability may reach a factor of a few. We 
have found only two ULX counterparts with si- 
multaneous X-ray — optical data. Below we de- 
scribe the observations and the SED model, and 
compare the observed and model SEDs. 

2. OBSERVATIONS AND DATA 
REDUCTION 

The catalog of ultra-luminous X-ray sources 
of Swartz et al. 2004 [29] provides data for 
more than 150 objects, however, less than 50 
objects are identified with optical sources [30]. 
There are considerably less thoroughly studied 
point sources for which the spectral energy dis- 
tributions are known in a wide range of wave- 
lengths ([21]). We have selected five ULXs that 



have reliable optical identifications with point 
sources: Holmbergll X-l [31], NGC6946 ULX- 

1 [19], NGC1313 X-l [32], NGC 1313 X-2 [33] 
and NGC 5408 X-l [34]. They are bright objects 
located in nearby galaxies (up to 5 Mpc) that 
have various spectral indices in the UV-optical 
range. 

Table 1 lists the basic information on the 
studied ultra-luminous X-ray sources, including 
the distances to the host galaxies, the adopted 
values of interstellar extinction, the X-ray lumi- 
nosities of the objects and the spectral indices in 
the UV-optical range. 

2.1. Optical data 

We use the HST archive data. The observa- 
tions were carried out with the Advanced Camera 
for Surveys (ACS), the Wide Field and Planetary 
Camera 2 (WFPC2) and the Wide Field Camera 
3 (WFC3) (the last one was used only for the 
observations of NGC 5408 X-l). 

For the objects NGC 1313 X-l, NGC 1313 X- 

2 and NGC 5408 X-l we used the photometric 
results from [21, 24, 32]. We selected quasi- 
simultaneous observations with a maximum set 
of filters. For NGC 1313 X-l we took the 
flux measurements in filters F330W, F435W, 
F555W and F814W taken on November 17, 2003, 
from [32]. For NGC 1313 X-2 the photometric 
measurements in filters F330W, F435W, F555W 
and F814W on November 22, 2003, were taken 
from [21]; for NGC 5408 X-l we used the re- 
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Table 1. Observational parameters of ULXs. The columns give: D — distance to the host galaxy in Mpc; Ay — total 
extinction, obtained from nebula spectroscopy; Ay G — extinction in the Galaxy; Lx — extinction-corrected X-ray 
luminosity in the 0.3-10 keV range in units of 10 ergs - ; a — UV-optical spectral index, F v <x v~ a 



Name 


D 


A v 


A 7 

A v G 


L x 


a 


Holmbergll X-l 


3.39 1 


0.23 


0.10 


6.0 8 


-1.07 ±0.05 


NGC6946 ULX-1 


4.5 2 


1.34 


1.14 


3.8 8 


-1.78 ±0.25 


NGC 1313 X-l 


4.13 3 




0.34 


15.7 9 


-0.81 ±0.03 


NGC 1313 X-2 


4.13 3 


0.40 5 


0.26 


5.4 9 


-1.34 ±0.04 


NGC 5408 X-l 


4.8 4 


0.25 6 


0.21 


10.2 10 


-1.32 ±0.07 



^arachentsev et al. 2002 [36], 2 N. A. Tikhonov [37], 3 Mendez et al. 2002 [38], 4 Karachentsev et al. 2002 [39], 
5 Grise et al. 2008 [40], 6 Kaaret & Corbel 2009 [41], 7 Schlegel et al. [42], 8 Swartz et al. 2011 [43], 9 average value 
based on the data from Feng & Kaaret 2006 [35], 10 Grise et al. 2012 [24] 



suits of the photometry in filters F225W, F336W, 
F547M, F845M, F105W and F160W on May 15, 
2010, from [24]. 

For Holmbergll X-l and NGC 6946 ULX-1, 
which were observed with HST in the far UV 
range, we performed photometry of the archive 
data. The information on the observation date 
and the used cameras and filters is listed in Ta- 
ble 2. 

For the observations obtained with the Ad- 
vanced Camera for Surveys (ACS), we selected 
the integral and distortion-corrected images in 
drz format. The standard initial reduction of the 
data from the ACS includes the procedures of 
bias and dark frame subtraction, flat-field cor- 
rection, accounting for "bad" columns and "hot" 
pixels, and cosmic-hit removal. The photometry 
of the WFPC2 data was performed on the images 
in the cOf format. The masking of defective pix- 
els, sky background subtraction and removal of 



cosmic hits were performed using the specialized 
program package HSTPHQT1 . 1. 

The fluxes of the optical counterparts of the 
X-ray sources were measured by means of aper- 
ture photometry using the APPHOT software 
package, which is based on IRAF. The aper- 
ture radii were selected in the / . / 15-0'. / 2 range. 
The main uncertainty in the point source flux 
measurements is due to the presence of very 
patchy nebulae. Following the paper by Kaaret 
et al. [19], we estimated the contribution of the 
nebula to the total flux in the aperture of the 
object and, based on this value, calculated the 
point-source flux error in every filter (Table 2). 

The aperture corrections for the data ob- 
tained with ACS/SBC were accounted for using 
the synphot program package. For WFPC2 and 
ACS/WFC the corrections were determined by 
photometric measurements of three to five single 
stars in the apertures varying from 0'.'15 to tf!5. 
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Table 2. HST data and results of photometry. The fluxes (given in units of 10 -17 erg cm~ 2 s _1 A -1 ) are 
corrected for interstellar extinction (Table 1). Total errors are given with allowance for the uncertainty associated 
with the subtraction of the contribution from the nebula. A p i vo t — the pivot wavelength of the filter in Angstroms 



Name 


Date 


Instrument / filter 


Apivot 


Flux 


NGC6946 ULX-1 


2008 May 1 


ACS/SBC/F140LP 


1533 


32.0±2.7 




2001 June 8 


WFPC2 /PCI /F450W 


4557 


2.3±0.6 




2001 June 8 


WFPC2/PC1/F555W 


5443 


1.09±0.27 




2001 June 8 


WFPC2/PC1/F814W 


7996 


0.27 ±0.05 


Holmbergll X-l 


2006 November 27 


ACS/SBC/F165LP 


1758 


33. Oil. 7 




2009 Febrary 9 


WFPC2 /WF2 /F336W 


3359 


4.62±0.18 




2007 October 3 


WFPC2/PC1/F450W 


4557 


1.98 ±0.18 




2007 October 5 


WFPC2/PC1/F450W 


4557 


1.88 ±0.17 




2007 October 9 


WFPC2/PC1/F450W 


4557 


1.96 ±0.18 




2007 October 3 


WFPC2/PC1/F555W 


5443 


1.06 ±0.07 




2007 October 5 


WFPC2/PC1/F555W 


5443 


1.07 ±0.07 




2007 October 9 


WFPC2/PC1/F555W 


5443 


1.07 ±0.07 




2006 December 30 


ACS/WFC/F814W 


8060 


0.301 ±0.012 



The so-called CTE effect may influence the re- 
sults of the measurements considerably, due to 
the fact that when reading the CCD, the effi- 
ciency of the charge packet transfer from pixel 
to pixel is not equal to 100%. When calculat- 
ing the CTE corrections, we used the internet 
resource CTE Tool ^l 1 , as well as the algorithm 
described in section 5.1.5 of the ACS Data Hand- 
book [44]. We have applied the CTE corrections, 
its value does not exceed 0™11. 

The flux within a given filter was corrected 
for interstellar extinction with the assumption 
that the spectra of the objects are described by 
the power law Fx oc A -3 . This procedure was 



1 http : //www. stsci . edu/hst/wfpc2/ sof tware/wfpc2_ 
cte_calc .html 



performed using the calcphot command in the 
synphot program package. 

The interstellar extinction for the objects 
NGC6946 ULX-1 and Holmbergll X-l was es- 
timated based on the Ha/H/3 flux ratio in 
the nebulae. For the photoionized plasma in 
widely varying conditions (at temperatures of 
5000 ± 20 000 K and electron densities of 10 2 ± 
10 6 cm" 3 ), Ha/H(3 is equal to 2.86 with an ac- 
curacy of about 5% [45] . In the case of collisional 
excitation this ratio may reach the value of 3.2, 
and then the ratio may be used for more 

reliable extinction estimates. 

We used the optical spectrum of the nebula 
surrounding Holmbergll X-l (taken on Febru- 
ary 28. 2011) from the SMOKA archive [46]. 



7 



To estimate the possible contribution of shocks 
to the observed Ha/H(3 ratio, we analyzed 
the lines of the nebula. The analysis of the 
[S II] AA6716,6734 doublet shows that the line 
widths of the doublet (corrected for the instru- 
mental profile 4.35 ± 0.15 A) amount to 92 ± 
8 km/s. This value agrees with the results of the 
velocity gradient measurements for Hell A4686 
in [17], which might indicate the presence of 
shock waves. However, the Ha/H/3 and H<y/H(3 
line ratios yield values of extinction that are con- 
sistent with each other and equal to Ay ~ 0.23. 
Therefore, the shocks contribution is fairly small. 
In our estimates of Ay we used the extinction 
curve of Cardelli et al. [47] with Ry = 3.1. The 
obtained value of extinction is in agreement with 
the results of [16]. 

For NGC 6946 ULX-1 we have performed sim- 
ilar extinction estimates based on the observed 
fluxes in the lines of the nebula reported in [16]. 
The Ha/H/3 ratio results in Ay ~ 1.55, whereas 
Hj/H/3 yields Ay ~ 1.34. Dunne et al. [48] have 
discovered broad emission line components with 
a velocity dispersion of about 250 km s _1 on the 
Echelle spectra of the object. It seems that in the 
case of the nebula around NGC 6946 ULX-1, the 
shocks affect the observed line ratios. For this 
object, we have adopted the value Ay ~ 1.34, 
measured from the H^f/H/3 line ratio. Extinc- 
tion corrected point-source fluxes are presented 
in Table 2. 

The values of interstellar extinction for 
NGC 1313 X-l, NGC 1313 X-2 and 



NGC 5408 X-l were taken from [40-42] (Ta- 
ble 1). In the case of NGC 1313 X-l, we used the 
Galactic value of extinction, since there are no 
extinction measurements based on the Balmer 
decrement of the nebula for this object. 

Figure 1 shows the spectral energy distribu- 
tions for the ultra-luminous X-ray sources. The 
lines show the results of the F v cx v~ a power- 
law approximation of the fluxes (Table 1). The 
fluxes are corrected for interstellar extinction and 
reduced to the Holmbergll X-l distance. 

For Holmbergll X-l, we used the data ob- 
tained at different epochs. The fluxes in filters 
F450W and F555W were averaged over three ob- 
serving dates. The best approximation result 
is the spectral index a = —1.07 ±0.05. In the 
case of NGC 6946 ULX-1 we used only optical 
data (a = -1.78 ± 0.25). The ultraviolet point 
falls out from the power law matching the op- 
tical data. For the remaining three sources we 
approximated the fluxes in all available filters. 
Obtained a are given in the Table 1. We used 
the 68% confidence interval in all of these esti- 
mates. 

2.2. X-ray data 

In this study we focus on the ULX spectra in 
the UV and optical ranges. Here, we take into 
account the X-ray data for illustrative purposes 
only. We therefore find it appropriate to use pub- 
lished results of the X-ray fits with simple mod- 
els. We selected the published models that best 
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Figure 1. Spectral energy distributions of the ULXs. The fluxes are corrected for interstellar extinction 
and reduced to distance D = 3.39 Mpc of the host galaxy of Holmbergll X-l. Simultaneous data are shown 
by same-type symbols. The spectra are shown from top to bottom: (a) NGC6946 ULX-1 (filled symbols), 
NGC5408 X-l, NGC 1313 X-2; (b) Holmbergll X-l, NGC 1313 X-l. With the exception of two objects 
(NGC6946 ULX-1 and Holmbergll X-l), we show only the simultaneous data to avoid overloading the 
picture. The lines depict the results of the power-law (F v oc v~ a ) approximation of the fluxes. 



match the observed spectra. All of these spec- 
tra have a sufficiently large total exposure time 
(> 10 ks). 

The Compton scattering plays an important 
role in the formation of the ULX X-ray spectra at 
energies greater than ~ 1 keV (e.g., [13]). We will 
model the X-ray data with the Compton scatter- 
ing in our next paper. Since we do not consider 
Comptonization here, we cut off the X-ray spec- 
tra above 1.5 keV. 

For Holmbergll X-l we selected the 
BB+DISKPBB model with p = 0.5 from [49], 
which was used to describe the observations 
from April 15, 2004. For NGC 6946 ULX-1, 
there are considerably less data available in 
the literature than for the other objects. We 
took the results of the DISKBB+POWERLAW 



model approximation of the September 7, 2001 
spectra of this source from [50]. For the two 
objects in the NGC 1313 galaxy we used the 
data from [35]. In the case of NGC 1313 X-l, 
we took the parameters of the models describing 
the X-ray spectra from October 17, 2000, 
January 17, 2004 and November 23, 2004. For 
NGC 1313 X-2 we used the results of the above 
authors for December 25, 2003 and August 23, 
2004. 

NGC 5408 X-l is the only object in our list 
for which three sets of simultaneous observations 
were obtained in 2010 in the optical (HST) and 
X-ray (Chandra) ranges. The results of these 
observations are published in the paper [24]. 
For our study, we took the data from May 15, 
2010, described by the DISKBB+POWERLAW 
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model. accretion rate; R m is the inner radius of the ac- 

In section 5 we separately compare the re- cretion disk; k is the Thomson opacity. Such an 

suits of the approximation by our model and by accretion regime is realized in the only known su- 

DISKIR for Holmbergll X-l. To this end, we peraccretor SS433, where Mq reaches the value 

re-analyzed the observations of Holmbergll X- ~ 10~ 4 MQ/year ~ 300MEdd [9, 51]. The su- 

1 obtained on April 15, 2004 with the XMM- percritical properties of the disk manifest them- 

Newton PN detector (OBSID 0200470101). The selves below the spherization radius: 

reduction was carried out using the SAS 12 pro- Mqk 

Rsp ~ -S— ■ ( 2 ) 

gram package. After event filtering with FLAG = 

PATTERNS 4, we selected the least noise- Within R sp , the disk is locally Eddington [26]: 
contaminated time interval of about 14 ks. Since the force of g ravit y is balanced by the radiation 
the object was situated close to the CCD gap, pressure at its every point. As a result, the disk 
we used the aperture with a radius of about 40" becomes geometrically thick (H/R ~ 1). Here- 
for extraction. The background was determined after > for certainty, we assume the half-opening 
in the ring around the object. The spectrum was an S le of the funnel of the accretion disk 6 f to be 
grouped in such a way that every bin contained ec i ual to 45 °> which corresponds to H/R = 1. 
at least 100 counts Because of the radiation pressure, the accret- 

ing matter begins to outflow from the surface of 
3 THE MODEL ^ e su P ercr itical disk in the form of wind, which 

leads to a gradual decrease of the accretion rate 
3.1. Shakura-Sunyaev supercritical disk with radius [26]: 

R 

The well-known paper of Shakura and Sun- M(R) = ^-M) (3) 

yaev [26] introduces the accretion disk model In tm ^ thig leadg to a decrease in the amoun t of 
(a-disk). The same paper considers the case of the released grav itational energy, and the result- 
supercritical accretion, which occurs when the ing digk mminosity becomes constrained by the 
rate of matter inflow into the disk is value 

L tot = L Edd (l + alnrho), (4) 



GMbh ck 



where 



where L Ed d ^ 1.3 x 10 39 Mi ergs 1 is the jiV D 

iviQ n sp 

m = — = t, — (5j 

Eddington luminosity for a black hole with M E dd ^in 

the mass Mbh expressed in units of 10 Mq; is the dimensionless initial accretion rate, and 

-^Edd = 2 x 10 19 Mio gs -1 is the corresponding a < 1 [8, 26, 52]. In this paper we adopt a = 1. 
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-^Edd) is released in the standard disk above R sp , 
Z/Eddlnmo is released in the supercritical disk. 

In the case of SS433, the relation (4) leads to 
the luminosity Ltot ~ 10 40 ergs -1 , which corre- 
sponds to the measurements in [53, 54] and co- 
incides with the independent estimate from [55]. 
The maximum of this luminosity lies in the ul- 
traviolet part of the spectrum. In the X-ray 
range, where the main contribution to the ob- 
served flux is given by the emission of the cooling 
relativistic jets, SS 433 is a relatively weak object 
(L x ~ 10 36 ergs" 1 ) [9]. 

The orientation of SS 433 is such that we ob- 
serve it close to the orbital plane and thus do not 
see the deep parts of the wind funnel of the super- 
critical disk. Initially, all of the energy release of 
SS 433 originates in the accretion disk in the form 
of hard radiation, which then thermalized in the 
powerful wind of the supercritical disk [9]. For 
an observer who sees the funnel of the supercrit- 
ical disk, SS 433 would have an X-ray luminosity 
of at least 10 40 ergs -1 . Ultra-luminous X-ray 
sources also demonstrate similar luminosities in 
the X-ray range [1, 29]. 

A further increase in the brightness of the su- 
percritical disk for the observer that can see the 
funnel is probably due to the geometrical colli- 
mation of radiation. In the simplest case, the 
radiation collimation factor is determined by the 
relation B = 2ir/Qf, where 0/ is the solid open- 
ing angle of the supercritical-disk funnel [7]. For 
the angle f ~ 45°, B ~ 3. 

In their paper [26] , Shakura and Sunyaev have 



pointed out the key role of the optically thick 
wind outflowing from the surface of the super- 
critical disk in the formation of UV and opti- 
cal spectra. This idea was developed further by 
many authors (e.g., [8]). 

We believe that the outflowing wind forms a 
funnel that has an approximately conical shape. 
Hereafter, we will call it "funnel" or "wind funnel". 
The existence of a funnel is due to the presence of 
an angular momentum both in the accreting and 
outflowing gas, which results in the formation of 
a region with a low density of matter near the 
symmetry axis of the system. The presence of a 
funnel both in the supercritical disk and in the 
wind outflowing from it is confirmed by direct 2D 
RHD simulations [27, 28]. 

The opening angle of the wind funnel is de- 
termined by the ratio of the Keplerian rotation 
velocity of the disk matter at a given radius to its 
outflow velocity, which is determined by the ra- 
diation pressure. The velocity of matter outflow 
from the surface of the supercritical disk should 
be close to virial [8, 26]; in this case we can ex- 
pect a wide funnel with the half-opening angle 
~ 30° -7- 60°. In this paper we adopt the open- 
ing angle of the wind funnel equal to the opening 
angle of the supercritical disk Of = 45°. 

Figure 2 shows the schematic model of the su- 
percritical disk with a wind funnel . Three main 
components are illustrated: the thin standard 
disk beyond the spherization radius, R > R sp , 
the supercritical disk (R m < R < R sp ), and the 
wind funnel, which spreads out to the radius of 
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the photosphere Rph/ sinOf. The photospheric 
radius confines the area of the optically thick 
wind; its size is determined by the geometry and 
velocity of the outflowing wind and by the depen- 
dence of the opacity coefficient on temperature. 
The exact calculation of i? pn is a rather difficult 
task; therefore, within the context of this paper, 
we consider it a model parameter. 

The typical sizes for a Schwarzschild black 
hole with a mass of 10 Mq and accretion 
rate ttiq = 300 are: R m = 3R g ~ 10 7 cm (R g 
is the gravitational radius), R sp ~3 x 10 9 cm 
(formula (2)). The radius of the pho- 
tosphere may be roughly estimated from 
the relation i? pn ~ M w k/4ttV w cos Of ~ 

2 x 10l2 (Ttra7i) _1 cm Q 56 ])> where m w ~ M 

is the mass loss rate, V w is the velocity of the 
wind at the level of the photosphere, which we 
assume to be 1000 kms -1 . 

3.2. Supercritical disk model (SCAD) and 
DISKIR model of Gierlinski et al. 

A supercritical disk can be observed as a lumi- 
nous X-ray source only if the angle between the 
line of sight and the axis of the disk is not greater 
than the funnel opening angle Of. In this case, 
an observer would see both the supercritical-disk 
radiation and the radiation of the wind funnel. 
In such an orientation, the central regions of the 
standard disk (R > R sp ), where the energy re- 
lease is maximal, may be hidden from the ob- 
server by the wind. We therefore assume that 



observer 




Tph 



TocR" 



3Rg Rsp 



Rph 



Figure 2. The model of a supercritical disk 
with a wind funnel. The figure shows the thin 
disk (R > -Rsp), the supercritical disk (R < R sp ), 
and the wind funnel constrained by the radius of 
the photosphere R p h / sin 8f . 

the radiation of the standard disk with the lu- 
minosity LecM heats the wind from below. This 
radiation is thermalized in the wind and thus in- 
creases its energy (see below). 

We assume that the radiation of the accretion 
disk and funnel is black-body-like. As was shown 
by Poutanen et al. [8], in a supercritical disk, the 
behavior of effective temperature with radius is 
determined by the relation 



T cxr 2 



(6) 



where r = R/R m . Such a temperature depen- 
dence leads to a power law spectrum F v cx u^ 1 . 

The temperature at the inner radius of the 
disk Tj n can be determined from the disk lumi- 
nosity 



L = 2 f 2nRaT A (R') dR' 

r' 

' sp 

= A,R m -aTj J ^_ . 



(7) 



Hereinafter, by R' = Rj sin Of we understand the 
distances measured along the disk and wind fun- 
nel surface; r' = R'/R m = r/sinOf. The 
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luminosity of the supercritical disk (R < R sp ), perature of the Comptonizing electrons), and 

the mass of the black hole, and the initial ac- L c /L& (ratio of the Compton tail luminosity to 

cretion rate are bound by the relation Lbol = the luminosity of the accretion disk), which de- 

a ^Edd In mo, a=l (4). Using (7) we find the tern- termine the shape of the spectrum in the inner 

perature regions of the disk. The parameters ri rr (radius 



-^Edd . a 
tt smt/f 

47raR in 2 . 



4 



I of the inner-disk area covered with "hot flow") 

(8) 

and f m (fraction of the Comptonized radiation 



A noticeable portion of the supercritical-disk that ends U P in the inner disk ) are responsible 

,. ,. , , , j , • j ri for the reprocessing of radiation in the r<r; rr re- 

radiation may be absorbed by the wind, tol- r e> 

, i i • • . i r t gion. With the exception of the L c /La ratio, the 

lowed by re-emission at lower frequencies. In ° ^ w u ' 

. , on at-, i i • • •, , ,i last 6 parameters have almost no effect on the 

this respect, the bCAD model is similar to the ^ 

model proposed by Gierlinski et al. (2009) [23], s i mu lated optical fluxes. 

DISKIR, which is currently used to explain the For the temperature T of the outer disk 

spectral energy distributions of ULXs [20, 21, 24]. heated h ? the radiation from the central parts, 

DISKIR is based on the standard "multi- Gierlinski et al. adopt Tocr"*. In the case of 

temperature disk" model (DISKBB), but takes heatin S of the standard (flat) accretion disk by 

into account the effects of Comptonization and the central source ' the temperature of the irra- 

irradiation of the disk by hard photons. The diated surface of the disk wil1 be determined by 

3 

j. i • ii ■ , , T oc r~4 the same as at gravitational energy re- 

mner disk is covered by a semi-transparent gas ' ° bJ 

cloud ("hot flow") up to the radius r irr . The lease ' Gie rlinski et al. suggested that the disk 

hot gas heats the outer parts of the accretion m ^ be war P ed ; in this case ' the temperature 

disk, which results in a considerable luminosity ^ de P end on radius as Tocr " 1 ' Since the tem " 

in the UV-optical range. Initially, this model has P e rature in the standard disk is determined by 

been applied to the low-mass X-ray source XTE the relation T K r_S > the reprocessed radiation 

J1817-330 [23] should dominate over the gravitational energy re- 

The DISKIR model contains 8 parame- leased in the disk at large radii, 
ters [23]. The observed UV, optical and 

near-infrared fluxes determine the parameters 3.3. The SCAD model parameters 
/out (fraction of reprocessed radiation) and r out 

(outer radius of the accretion disk). The X-ray In our model we assume that the relation 

part of the spectrum constrains the other six pa- Tctr~^ also applies for the wind funnel. Such 

rameters, including kT[ n , T (power law photon a dependence may be caused by the geometry of 

index for the Comptonized radiation) , kT e (tern- the wind (warped wind) or the presence of scat- 
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tering matter filling the funnel. 

Similar to [23], we determine the re-emitted 
flux by the following equation: 

/out-^tot / n \ 

where / out is the fraction of the bolometric flux 
reprocessed in the wind at the radius R', and L tot 
is the bolometric luminosity (4). The resulting 
spectrum of the wind funnel, as well as the spec- 
trum of the supercritical disk, will be power law, 
F v (x I/ -1 . 

The parameter / out is a combination of many 
geometrical and physical characteristics of the 
wind, such as the fraction of the accretion-disk 
radiation absorbed by the wind, the albedo, and 
the absorbed-radiation thermalization efficiency. 
The supercritical disk with luminosity -^Edd hi fho 
heats the wind from the side of the funnel, and 
the standard disk at R > R sp with luminosity 
Lgdd heats the same wind from below. 

As a result, the total flux of the supercritical 
disk with a wind funnel will be described by 

(r'sm9 f )~ 2 , 1 < r' < r sp ' 
F to t = vT\A ^-(l + lnmoV" 2 , 
k r sp ' < r' < r ph ' 
Thus, our model has seven parameters. The 
four main parameters that fully determine the 
shape of the model spectra are: black hole mass 
Mbh (-^Edd oq Mbh), initial normalized accretion 
rate tHq, fraction of bolometric flux reprocessed 
in the wind / out , and wind photosphere radius 
R' ph . Next two parameters are the opening angles 
of the supercritical disk and wind funnel, and in 




Wavelength, A 



Figure 3. Model spectral energy distributions 
(SEDs) for the supercritical disk with a wind 
funnel. Different grey tones correspond to 
different black hole masses Mbh for the accretion 
rates m = 30 (dots lines), 100 (dashed lines) and 
500 (solid lines) (a), and different values of / out 
for the photospheric radii R' ph /Ri n — 0.5 x 10 5 
(dots lines), 1.5 x 10 5 (dashed lines) and 3 x 10 5 
(solid lines) (b). The values of the fixed 
parameters are: / ou t = 0.1 and R' ph = 10 12 cm 
(a), Mbh = 15 M Q and rh = 200 (6). 

this paper, we assume them to be equal, 9f = 
45° . We set the seventh parameter a (4) equal to 
one, assuming within this paper that advection 
is absent in the accretion disk (see section 5). 

Figure 3 shows the model spectral energy dis- 
tributions. The fluxes were computed for the dis- 
tance D = 3 Mpc. Two flat areas on the spectra 
correspond to the power law F v cx v~ l . The radi- 
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ation of the supercritical disk dominates in the X- 
ray range. In the UV/optical range the observed 
spectra are mainly determined by the radiation 
of the irradiated wind. In the long-wavelength 
part, the flux is determined by the Rayleigh- 
Jeans spectrum of the wind funnel, F u oc v 2 . 

Figure 3a demonstrates the change in the 
spectra as a function of black hole mass Mbh and 
accretion rate mo- With increasing mass, the Ed- 
dington limit Lgdci) and therefore, the bolomet- 
ric luminosity of the accretion disk proportion- 
ally increase. The luminosity of the wind funnel 
also linearly grows with Mbh when / out = const. 
Therefore, on the whole, the radiation spectrum 
shifts to the domain of higher fluxes. 

An increase in the accretion rate at constant 
black hole mass leads to the same change in the 
geometric size of the supercritical disk and a de- 
crease in the gas and radiation temperature in its 
outer regions. For the temperature at the spher- 

ization radius, the following relation is valid: 

_ i 

T sp oc m 2 . The temperature decrease in the 
outer parts of the disk is due to the increase of 
the radiating area (i? sp 2 ocmQ). As a result, the 
flat part of the X-ray spectrum extends to the 
long- wavelength region (Fig. 3a). 

The dependence of the spectra on / out and 
R' ph at constant values of Mbh and mo is shown 
in Fig. 3b. The parameter / ou t mainly acts as a 
scale factor for the UV and optical parts of the 
spectrum, since the flux linearly grows with the 
fraction of absorbed and re-emitted energy. At 
the same time, the dependence of radiation tem- 



perature on the amount of energy reprocessed 
by the wind is fairly weak and proportional to 
/out 3 - For this reason, the wind temperature de- 
creases with decreasing / out , and the flat part 
of the UV-optical spectrum shifts slightly to the 
long- wavelength region (Fig. 3b). 

The radius of the photosphere determines the 
outer boundary of the blackbody emitting wind. 
The increase of the wind-funnel size at a con- 
stant irradiating flux leads to the decrease in the 
intensity and temperature of the wind radiation. 
Thus, the increase of the photospheric radius of 
the wind funnel extends the flat part of the UV 
spectrum into the optical range. 

4. RESULTS 

We have fitted the ULX SEDs with the 
SCAD model. As X-ray data, we used the pub- 
lished models (including absorption models) that 
best describe the observed spectra (section 2.2). 
Based on the parameters of these models, we 
have constructed the absorbed spectra. We used 
models instead of the observed spectra because 
the data approximation was carried out without 
XSPEC. 

When approximating the absorbed X-ray 
spectra with our model, we determined the new 
value of Nff. We used the extinction model 
from [57]. In the optical range we used the 
extinction-corrected fluxes. Unlike in the X-ray 
range, the optical interstellar extinction is usu- 
ally well-known from the spectroscopy of the neb- 
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ulae surrounding the objects. 

Using the relation Ay = 4.5 x W~ 22 N H [58], 
we have compared the extinctions in the optical 
and X-ray ranges (Tables 1 and 3). In all cases, 
the value of Ay obtained from the optical data 
is smaller than the Ay obtained from the X-ray 
spectra. The difference (-4y) X - r ay ~~ (^)opt is 
~ a"2 4- 0P8, the average value is ~ (P5. 

Table 3 shows the results of spectral energy 
distribution simulations for all of the five objects. 
The table includes the values of column density 
taken from literature (iV#) lit and the ones ob- 
tained by us, Nh, inner radius temperatures, 
kT[ n , black hole masses Mbh, normalized accre- 
tion rates mo, fractions / out of the bolometric 
flux of the accretion disk reprocessed in the wind, 
and radii of the wind- funnel photosphere R' ph . 
We do not quote the parameter errors, because in 
our approximations we used X-ray spectra, which 
are not direct observational data. The obtained 
values of Mbh correspond to stellar-mass black 
holes, and riiQ are comparable to the accretion 
rate in the SS433 system. The average values 
of parameters / out and R' ph for the five studied 
objects are 2.8 x 1CT 2 and 2.3 x 10 12 cm, respec- 
tively. 

Figures 4 and 5a show the observed spectral 
energy distributions and their approximation by 
the SCAD model. The optical fluxes corrected 
for interstellar extinction are shown with the er- 
ror bars. The black crosses and grey dots indicate 
the absorbed X-ray spectra reconstructed from 
literature. To avoid overloading the images, the 



spectra demonstrating object variability in the 
X-ray range are shown only for NGC 1313 X-l 
and NGC 1313 X-2 (grey dots). 

The curves show the model and its compo- 
nents in different spectral regions. The spec- 
tra forming in the supercritical accretion disk 
(dashed curve) and wind funnel (dot-dashed 
curve) together give the integral spectrum, 
shown by upper curve. At higher energies the 
solid line shows our model approximation of the 
X-ray data. The lower panels show the discrep- 
ancy between the observations (O) and the model 
(C). 

NGC 1313 X-l, NGC 1313 X-2. For 
NGC 1313 X-l and NGC 1313 X-2 the model 
SEDs describe both the optical and X-ray data 
well (Fig. 4a,b). The deviations do not exceed 
10-15%. A distinctive feature of NGC 1313 X-l 
is the very luminous (compared to the optical 
fluxes) X-ray spectrum, which is unusual for 
our sample. The value of / out for this object is 
more than one order of magnitude less than the 
average value in other objects. 

Holmberg II X-l . The model describes the en- 
ergy distribution of Holmberg II X-l fairly well 
(Fig. 5a). The only outlier is the flux in fil- 
ter F336W (Api vo t = 3359 A) that approximately 
corresponds to filter U, which may be due to a 
Balmer jump or object variability. Note that 
non-simultaneous optical observations were used 
in the approximation of this ULX. 

NGC 6946 ULX-1. NGC 6946 ULX-1 
(Fig. 4c) demonstrates a very low X-ray flux 
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Table 3. Model parameters of ULXs. {Nn) lit — column density taken from literature, Nh — column density 
obtained from our fit to the X-ray spectra; kTi n — inner disk temperature; Mbh — black hole mass; rho — initial 
normalized accretion rate; / ou t — fraction of the bolometric flux reprocessed in the wind funnel; _Rp h — radius of the 

photosphere 
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kT- ln 
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font 


R ' P h 




(x 10 21 cm' 




(x 10 21 cm- 2 ) 


(keV) 


(Mq) 




(x 10- 2 ) 


(x 10 12 cm) 


Holmberg II X- 1 


0.87 1 




0.94 


0.81 


17 


270 


4.1 


2.3 


NGC6946 ULX-1 


6.12 2 




4.3 


0.77 


20 


180 


4.2 


3.5 


NGC 1313 X-l 


2.5 3 




2.5 


0.82 


16 


150 


0.18 


1.6 


NGC 1313X-2 


2.7 3 




2.5 


0.99 


8 


200 


1.8 


1.3 


NGC 5408 X-l 


1.27 4 




0.97 


0.85 


14 


250 


3.8 


2.6 



^ajava et al. 2012 [49], 2 Berghea et al. 2008 [50], 3 Feng & Kaaret 2006 [35], 4 Grise et al. 2012 [24]. 



(both in the soft and hard ranges) compared to 
the UV and optical fluxes, which distinguishes 
this object from the other ULXs in our sample. 
To reconcile the model with the X-ray data, we 
excluded the areas with R < 5R m . We sug- 
gest that this object might be visible at an angle 
slightly larger than the opening angle of the wind 
funnel 6j. In this case, the innermost parts of the 
supercritical-disk funnel may be hidden from the 
observer by the outer parts of the wind. 

NGC5408 X-l. For NGC 5408 X-l (Fig. 4d), 
we were unable to achieve a good agreement be- 
tween the model and the X-ray data at energies 
higher than ~ 0.8 keV (~ 16 A), where the sim- 
ulated fluxes are twice as high as the observed 
ones. In our next paper we intend to include 
the Comptonization in the inner parts of the su- 
percritical accretion disk. Hot wind may cover 
the inner parts of the disk, absorbing and Comp- 
tonizing their radiation. We suppose that in the 



case of NGC 5408 X-l, this wind captures sub- 
stantially more disk radiation than in the other 
objects. We attribute the discrepancy between 
the model and the X-ray data to the presence 
of such a wind. In the optical range, our model 
reconstructs the observed fluxes fairly well. A 
slight excess is observed only in the IR range in 
filter F160W (Ap ivot = 15369 A). 

Thus, on the whole, the SCAD model gives 
good agreement between the observed and simu- 
lated energy distributions of ULXs in the wave- 
length interval from soft X-ray to optical. 

Figure 5 shows the SED approximation by the 
SCAD model (a) and by DISKIR (b) for Holm- 
berg II X-l. To correctly compare the approxi- 
mation results of the two models, we used the 
same spectra from April 15, 2004. In Fig. 5a, 
this X-ray spectrum is shown as the best model 
from [49], and Fig. 5b shows the spectrum pro- 
cessed by us (section 2.2). When using DISKIR, 
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Figure 4. The observed energy distributions of ULXs and their fit with the model of a supercritical disk 
with a wind funnel for NGC 1313 X-l (a), NGC 1313 X-2 (b), NGC6946 ULX-1 (c) and NGC5408 X-l (d). 

The grey dashed curve shows the spectrum formed in the supercritical accretion disk, the grey dash-dotted 
curve shows the spectrum formed in the wind funnel, and the grey solid curve is the integral spectrum. The 
black crosses (and grey dots in the cases of NGC 1313 X-l and NGC 1313 X-2) represent the absorbed X-ray 
spectra, and the black curve shows our model approximation of the X-ray spectrum. The lower panels show 
the relative deviations of the observed energy distributions from the simulated ones. 



the extinction in the X-ray region was recon- 
structed using the WABS model, and in the 
optical range, we approximated the extinction- 
corrected fluxes. 

Following [20, 23], we fixed the parameters 
of electron temperature, kT e = 100 keV, inner- 
disk radius, r m =1.1, and fraction of Comp- 
tonized radiation thermalized in the r < r- lxr 
region, /; n = 0.1. The remaining parame- 



ters were considered free. The best agree- 
ment between the model and observations 
with x 2 /dof = 340.0/299 was obtained with the 
values 

N H = (1.06 ± 0.16) x 10 21 cm" 2 , 

kT m = 0.12 ± 0.03 keV, T = 2.53 ± 0.03, 

L c /L d = 2.7 ± 1.7, / ou t = (3.3 ± 0.5) x 10" 2 , 
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Figure 5. Comparison of the SCAD (a) and 
DISKIR (b) models for Holmbergll X-l. 
Designations as in Fig. 4. 



log (r out ) = 3.23 ± 0.12 (in R in radii), 

norm = (7 ± 3) x 10 2 . 

Both models describe the X-ray and optical 
data equally well. The most significant differ- 
ence between the model energy distributions is 
in the slope of the X-ray spectrum in the ~ 0.3- 
1 keV region (12-40 A), where the model of the 
supercritical disk predicts lower (compared to 
DISKIR) fluxes and a flat spectrum. Both mod- 
els result in similar values of Nh, however, the 
value of this parameter in DISKIR is larger by 
~10%. This is caused by the necessity to com- 



pensate for the higher (compared to the SCAD 
model) fluxes in the soft X-ray range. 

Nevertheless, the models predict vastly dif- 
ferent temperatures at the inner disk radius. 
For Holmbergll X-l, the SCAD model gives 
kT- m ~ 0.8 keV, whereas DISKIR yields an es- 
timate of kT m ~ 0.1 keV. In turn, the low 
temperature at the inner radius leads to a 
black hole mass increase. Using normaliza- 
tion in the DISKIR model and a mean cosine 
of the angle between the disk and the line of 
sight < cos# >= 0.5, we derive R{ n ~ 10 9 cm, 
and the corresponding mass Mbh ~ 1000 Mq 
(L/L Edd ~ 0.05). The SCAD model gives 
M BH ^ 17 M for Holmbergll X-l. 

5. DISCUSSION 

A comparison of our model and DISKIR in 
the case of Holmbergll X-l (Fig. 5) shows that 
both models well describe the observed optical 
to soft X-ray spectrum. The DISKIR model has 
already been applied to Holmbergll X-l in the 
study [21]. Despite the fact that different obser- 
vational data was used by these authors and us 
(Fig. 5b), the parameters of the best SED model 
in both cases turned out approximately the same. 
In particular, for Holmbergll X-l, [21] derived 
the masses 400 and 600 Mq from the obtained 
inner-radius temperature kT m = 0.09-0.24 keV. 

For the other ULXs to which the DISKIR 
model has been applied, NGC 6946 ULX-1 [25] 
and NGC 5408 X-l [24], the authors obtained 
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similar kT in temperatures, 0.12 and 0.13 keV, re- 
spectively. Obviously, all of these temperatures 
correspond to the IMBH mass range; the [23] 
model gives large black hole masses for the ULXs. 
Indeed, the DISKIR model presumes a standard 
accretion disk, where the huge X-ray luminosi- 
ties of ULXs can be sustained only in the case 
of IMBHs. However, the recent data [5] reliably 
shows that ULXs are stellar mass black holes, 
the progenitors of which are massive stars with 
masses of 50-100 M . Our SED model, which 
is based on the SCAD approximation of Shakura 
and Sunyaev, presumes black hole masses of 10- 
20 Mq, which can well be achieved in the pro- 
cess of stellar evolution. Therefore, this model is 
compatible with the new data on the nature of 
ULXs. 

The shape of the spectrum in the soft X-ray 
region may probably serve as a critical test for 
distinguishing between the SCAD and DISKIR 
models. In this region, 12-40 A (~ 0.3 - 1 keV), 
the SCAD model demonstrates a flat spectrum 
{vF v ), whereas in the DISKIR model there can- 
not be any flat parts anywhere, as it is based on 
the standard a-disk. In DISKIR, when L c /L^ 
changes from 0.3 to 7, the maximum in the model 
spectrum shifts from 0.3 keV energies to 1.0 keV. 

In Figure 5, the differences between the 
SCAD and DISKIR models are not very con- 
spicuous, both the best models appear the same. 
This is due to the fact that the other param- 
eters of the models, such as Njj and disk lu- 
minosity, when combined in a certain way, can 



smooth out the differences between the models. 
Nevertheless, the absence of a flat spectrum in 
the DISKIR model and the presence of such in 
SCAD are essential. High signal-to-noise spectra 
may probably be used to differentiate between 
the two models. 

In this paper we assumed some of the model 
parameters to be constant, e.g., the opening an- 
gle of the supercritical-disk funnel Of and the 
opening angle of the wind funnel, which is formed 
by the wind of the supercritical disk, were set 
equal to 45°. Within the framework of our sim- 
ple approach, a change in the opening angle of 
the funnel will change the observed luminosity, 
and therefore, as a result, will somewhat re- 
determine the black hole mass, Mbh °c sinfly -2 . 
This equally applies to the parameter a (4), 
which we set equal to one, as would be the case 
in the absence of any advection in the accre- 
tion disk [26]. Here M B h oc (1 +alnm )~ 1 . At 
rriQ = 300, a decrease of a from 1 to the expected 
value ~ 0.5 ([59, 60]) will lead to a 1.75 times 
increase of the mass estimate. In this paper we 
also disregard the collimation of radiation by the 
funnel in the supercritical disk and wind, which 
can amount to ~ 2-3 if the funnel opening angle 
is 45-55°. Taking into account the collimation 
of the radiation toward the observer will lead to 
a decrease in the black hole mass estimate by 
about ~2 times. 

In the current version of the model we disre- 
gard Comptonization— a primary component in 
the formation of X-ray spectra. We assume that 
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the inner parts of the supercritical disk are cov- 
ered by a semi-transparent hot gas. By analogy 
with SS433, one can claim that such gas should 
also exist in the funnel. This gas may not become 
a part of the jets during their formation and col- 
limation, or, it is this gas that constitutes the 
relativistic jets of SS433 at the moment of their 
formation. The density of this gas may roughly 
vary according to the relation n oc r~ 2 ; this gas 
should leave the funnel. The presence of such gas 
will lead to the Comptonization of the radiation 
of the supercritical-disk inner parts, which are 
covered by this gas. 

Medvedev and Fabrika [61] have discovered 
an X-ray emission component in the XMM data 
of SS433 that has a flat spectrum. The au- 
thors suggested that this is Comptonized radia- 
tion, formed in the inner parts of the funnel and 
reflected by the funnel outer walls. An analogy 
with SS 433 gives us a reason to assume that the 
inner regions of supercritical disks may contain 
hot outflowing gas. After including Comptoniza- 
tion, which we plan to do in the next paper, we 
will be able to approximate the ULX spectra in 
the entire standard X-ray range. 

In this paper we assumed that the UV and 
optical radiation in ULXs is formed in the ac- 
cretion disks (standard, in the case of IMBHs, 
or supercritical in the case of stellar-mass black 
holes). However, this radiation may be formed 
at the surface of the donor. The optical spec- 
tra of ULX counterparts show that all the ob- 
jects ever observed have photosphere tempera- 



tures T > 30 000 K, and it is unlikely that we 
observe the donor's own radiation in the UV and 
optical ranges [10]. 

Let us assume that the IMBH heats the 
donor's surface and the latter overfills its critical 
Roche lobe. The size of the donor in units 
of orbital separation [62] is r^/a ~ O^g 1 / 3 , 
where q = Md/Mbh and Md the donor mass. 
If we adopt the X-ray isotropical luminos- 
ity of ULXs L x = 10 40 L4o erg/s, and the 
donor surface albedo A = 0.5, we derive 
the bolometric luminosity of the heated donor 
L bol ~ AL x {r D /a) 2 /A ~ 1.2 x lO^L^oii erg/s, 
where <?o.oi the mass ratio in units 0.01. This 
luminosity corresponds to the bolometric 
magnitude of the donor Mboi ~ —4.0. For a 
star with the temperature, T > 30 000 K the 
bolometric correction is BC < —3.0. We find 
the absolute magnitude of the heated donor 
My > —1. This is substantially smaller than 
the counterpart luminosity observed in ULXs 
(My ~ — 6±1 mag [21]). Therefore, the scenario 
of the formation of optical and UV spectra on 
the irradiated donor surface must be excluded. 

6. CONCLUSIONS 

In this paper we have described a model of en- 
ergy distribution in supercritical accretions disks, 
which is based on the Shakura-Sunyaev (1973) 
supercritical disk conception. We applied this 
model to ultra-luminous X-ray sources. In this 
model, the disk becomes thick at distances from 
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the center less than the disk spherization ra- 
dius, H/R ~ 1; the temperature depends on ra- 
dius as T oc r -1 / 2 . Here, the disk luminosity is 
-^bol — -^Eddhimo, and also, in this region, the 
powerful wind that forms a funnel above the disk 
emerges. At distances greater than the spheriza- 
tion radius, the disk is thin; its total luminosity 
is equal to the Eddington luminosity LecM- 

We believe that for an observer that can see 
the entire funnel of the supercritical disk (who is 
situated close to the accretion-disk axis, < 9f), 
the thin accretion disk at the radii r > r sp is 
invisible, as it is hidden by the extended wind of 
the supercritical disk. The thin disk irradiates 
this wind from below. From the side of funnel, 
the wind is irradiated by the supercritical disk. 
In this paper we adopt the following dependence 
of the irradiated wind temperature on distance: 
T oc r -1 / 2 ; this dependence may be due to both 
the curvature of the wind, and the presence of gas 
in the funnel, in which jets possibly spread. Both 
these effects cause an increase in the efficiency of 
irradiation of the wind funnel walls. 

In this paper we do not consider the Comp- 
tonization of the radiation emerging from the 
central regions of the supercritical disk funnel. 
Accounting for Comptonization is absolutely nec- 
essary in this model, but we restrict this study 
to the analysis of the spectra from ~ 1.5 keV to 
the lower energies. By analogy with SS433, we 
assume that there should exist a hot wind in the 
inner parts of the funnel; this is the gas of the 
forming jets (and also the gas that contributes 



to the formation of the jets, but does not be- 
come part of them) . This wind covers the inner 
parts of the supercritical disk and Comptonizes 
the radiation from these parts. 

The described supercritical disk model is 
technically similar to the DISKIR model [23], 
which was initially developed for LMXBs, but 
has already been used for the analysis of ULX 
SEDs. The models differ in disk type (standard- 
supercritical) and in the fact that in one model, 
the disk is irradiated at the distances r > rj rr , 
and in the other, the wind is irradiated at r > r sp . 
We believe that distinguishing between the mod- 
els is possible in the soft X-ray range (~ 0.3 — 
1 keV), since the SCAD model has a flat region in 
the spectrum, whereas DISKIR cannot have flat 
regions anywhere, as it is based on the standard 
a-disk. Since the differences between the models 
can be masked by varying the parameters (e.g., 
Njj), obtaining X-ray spectra with high signal- 
to-noise ratios is a necessary condition for distin- 
guishing between the models. 

We see a rather important difference be- 
tween the models in the fact that DISKIR, when 
applied to ULXs, explains their energy distri- 
butions only in the case of black holes with 
masses of several hundred solar masses (IMBHs), 
whereas the SCAD model offers fairly adequate 
black hole masses ~ 10 Mq, which can be ob- 
tained in the course of evolution of massive stars. 
Recent results of ULX studies favor the nature of 
these objects as supercritical accretion disks with 
stellar-mass black holes, but not IMBHs. 
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To adequately describe the ULX SEDs in any 
model, in particular, in SCAD, one needs to ob- 
tain simultaneous observational data in the X- 
ray to optical range, and most importantly, in 
the soft X-ray and UV spectral regions. Since 
in both the SCAD and DISKIR models the UV 
and optical spectra are formed as a result of ir- 
radiation of the outer regions by the radiation 
from the inner regions, the observations must be 
simultaneous in all of these ranges. 
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